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ABSTRACT 

Core-collapse supernovae (SNe) are the spectacular finale to massive stellar evolution. In this 
Letter, we identify a progenitor for the nearby core-collapse SN 2012aw in both ground based near- 
infrared, and space based optical pre-explosion imaging. The SN itself appears to be a normal Type 
II Plateau event, reaching a bolometric luminosity of 10 42 erg s _1 and photospheric velocities of 
^11,000 km s" 1 from the position of the H/3 P-Cygni minimum in the early SN spectra. We use an 
adaptive optics image to show that the SN is coincident to within 27 mas with a faint, red source 
in pre-explosion HST+WFPC2, VLT+ISAAC and NTT+SOFI images. The source has magnitudes 
i^555IF=26. 70+0.06, F814JE=23.39±0.02, J=21.1±0.2, A=19. 1+0.4, which when compared to a 
grid of stellar models best matches a red supergiant. Interestingly, the spectral energy distribution 
of the progenitor also implies an extinction of Ay >1.2 mag, whereas the SN itself does not appear 
to be significantly extinguished. We interpret this as evidence for the destruction of dust in the SN 
explosion. The progenitor candidate has a luminosity between 5.0 and 5.6 log L/Lq, corresponding 
to a ZAMS mass between 14 and 26 M Q (depending on Ay), which would make this one of the most 
massive progenitors found for a core-collapse SN to date. 

Subject headings: supergiants — supernovae: general — supernovae: individual (SN2012aw) — galax- 
ies: individual (M95) 



1. INTRODUCTION 

That massive red supergiants explode as Type 
IIP (Plateau) supernovae (SNe) at the end of 
their lives has been clearly estab l ished (for example. 
Van Dvk et al.l 120031 ISmartt et all 12004 iMattila et all 
20081 iFraser et al.l 120111) , and confirmed by the disap- 
pearance of their progenitors (jMaund fc Smarttl 120091 ) . 
What remains of great interest, however, is understand- 
ing how SN progenitor properties (eg. mass, radius, 
metallicity) correlate with the explosion characteristics 
(kinetic energy, synthesised 56 Ni, luminosity etc). The 
continued identification of progenitors of nearby SNe, 
coupled with follow-up observations of the SNe, can help 
shed light on this link. In this Letter we present an anal- 
ysis of the progenitor of SN 2012aw, which is the coolest, 
and probably the most massive progenitor of a Type IIP 
SN found thus far. 

SN 2012aw in M95 was dis covered independen tly by 
several amateur astronomers ()Fagotti et al.l 120121) . with 
a first detection on 2012 Mar 16. 9. A non-detection at a 
limiting mag nitude of R > 20.7 (|Poznanski et all 120121) 
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on Mar 15.3 sets a rigorous constraint on the age of the 
SN at discovery of < 1.6 days, an d we adopt an explosion 
epoch of Mar 16.0 UT (± 0.8 d). IMunaril (|20T1 obtained 
a spectrum of the SN on March 17.8, which sh owed a fea- 
turel e ss blue continuum; subsequent spectra (jltoh et al.l 
120121: ISiviero et aLll2012[ ) suggested the SN was a young 
Type IIP, although this classification is yet to be con- 
firmed by the emergence of a plateau in the lightcurve 
at 3-4 weeks after explosion. A candidate progenitor for 
SN 2012aw w as identified in archival Hubble Space Tele- 
scope data by [Elias-Ros a et al.l (|2012D and IFraser et al. I 
(2012), who both suggested that a faint red supergiant 
was the likely precursor star. 

The host galaxy of SN 2012aw, M95, is a barred 
and ringed spiral gal axy with an inclination of 55°. 
iFreedman et al . (2001) measured a Cepheid-based dis- 
tance to M95 of 10 Mpc (corresponding to p. = 30.0±0.09 
mag), which we adopt in all of the the following. This dis- 
tance is consistent with that obtained from the Tip of the 
Red Giant Branch (Ri zzi et all 12007). We also adopt a 
foregroun d (Milky Way) extinc tion of E(B -V) = 0.028 
mag from Schlcgel et al.l (|1998D . Usi ng the line strengths 
for the closest Hn region reported in lMcCall et al.l (|1985D 
we calculate a metallicity of 12 +log[0/Hl = 8.8+0.1 us - 
ing the [Om]/[Nii] relation of iPettini fc Pagell (|2004D . 
From the rec ent study of the radial metallicity gradi- 
ent in M95 bv lPilvugin. Thuan fc Vilchezl ([20061 ). at the 
position of SN 2012aw we estimate a metallicity of 12+ 
log [O/H] = 8.6±0.2. Hence we suggest that the metal- 
licity of the progenitor of SN 2012aw is approximately 
solar (for comparison, Milky Way Hn regio ns have a typi- 
cal m etallicity of 12+ log [O/H] = 8.7+0.3; iHunter et al.l 
I2009D , albeit with the caveat that metallicity calibrations 
ha ve significant uncert ainties (as discussed in this context 
bv ISmartt et al.l 120091 ) . 
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SN 2012aw wa s detected in X-rays wit h the Swift X- 
ray Telescope bv llmmler fe Brownl (|2012j), at a luminos- 
ity of L x = 9.2 ± 2.5 x 10 38 erg s" 1 . This is at the 
upper end of the range of Type IIP X-ray lu minosities 
(0.16 - 4 x 10 38 erg s" 1 ; iChevalier et al.|[200l ). although 
still much lower than the typica l values seen fo r Type 
Iln SNe (0.2 - 1.6 x 10 41 erg s'^ lFox et aT1[2000l ) which 
are interacting with significant amounts of circumstellar 
material . SN 2012aw was also detected in radio obser- 
vations (jYadav et al.ll2012t iStockdale et al.ll2012fl with a 
flux at 22 GHz of 2xl0 25 erg s" 1 Hz" 1 at +1 week after 
explosion, rising to 4xl0 25 erg s Hz -1 at +2 weeks. 
These values are simi lar to those of the sam ple of Type 
IIP SNe presented bv lChevalier et al.l (|2006ft . 

2. SUPERNOVA CHARACTERIZATION AND FOLLOWUP 

Our collaboration^] commenced an intensive spectro- 
scopic and photometric follow-up campaign for SN 
2012aw immediately after the discovery. Full coverage 
of the SN will be published in a future paper, here we 
present a limited set of optical observations from first 
two weeks after explosion to characterize the initial SN 
evolution. 

Photometry was obtained for SN 2012aw with the 
Asiago Observatory Schmidt telescope in the Landolt 
BVRI system. The data were pipeline reduced, and 
PSF-fitting photometry was performed on th e images us- 
ing the QUBA pipeline (jValenti et al.l 120111 ). A pseudo- 
bolometric (ie. BVRI) lightcurve was constructed by 
integrating the observed SN flux over the optical fil- 
ters. The resulting lightcurve is shown in Fig. [TJ to- 
gether with those of three other Type IIP S Ne for com- 
parison (SN 1999em , lElmhamdi etail [20031: SN 1999m , 
iLeonard et aT1[200l SN 2005cs. iPastorello et"all[200l . 
The luminosity of SN 2012aw is comparable to that of the 
SN 1999em, placing it firmly in the continuum of normal 
Type IIP SNe, as opposed to sub-luminous events such 
as SN 2005cs. 

Spectra of SN 2012aw were obtained from the Asi- 
ago 122cm telescope + Boller&Chivens Spectrograph + 
300tr/mm, and the Nordic Optical Telescope + ALFOSC 
+ Gr#4. Spectra were reduced, extracted, and flux and 
wavelength calibrated with the QUBA pipeline. The se- 
quence of spectra from the first two weeks after explosion 
is presented in Fig. [TJ Without any adjustment for ei- 
ther Milky Way or host galaxy extinction, the slope of 
the SN 2012aw spectrum (epoch 3.9 d) is almost iden- 
tical to that of SN 1999em at +4 d (jBaron et al.l 120001: 
IHamuv et al1l200l . H aron et all found that a synthetic 
spectrum fit with T-qb = 11000 — 13000 K and an extinc- 
tion of E(B — V) ~ 0.05 mag provides a good match 
to the continuum slope and Hi lines. A similar result 
would necessarily be determined for SN 2012aw given 
its striking similarity to SN 1999em. The Ha absorp- 
tion minimum is at a velocity of v p hot ~ 13000 km s _1 , 
again almost identical to the minimum for SNe 1999em 
and 1999gi. Fig. [TJ shows the early spectra of these three 
Type IIP SNe illustrating the similarity, also shown is the 
low energy, faint Type IIP SN 2005cs (jPastorello et al.l 
12009ft which has a Ha minimum at a much lower velocity 
of 6500 km s _1 . This clearly shows that SN 2012aw is 

9 http:/ /graspa.oapd. inaf.it/index. php?option=com_content&;vi 
ew=article&id=68&Itemid=93 



not a low energy explosion. 

A high resolution spectrum of SN 2012aw was obtained 
with the Telescopio Nazionale Galileo + Sarg on Mar 29. 
The two components of the Nai D doublet were detected 
at 5907.998, 5914.004 A , which is consistent with the 
recessional velocity of M95. The equivalent width (EW) 
of the two components were 0.286 a nd 0.240 A respec- 
tively. From the empirical relation of IMunari fe Zwitterl 
((19971 ). we estimate E(B-V) = 0.10±0.05 mag for SN 
2012aw. We note however, that even if dust has been 
photo-evaporated in the explosion (see Sect. [4j, the Nai 
absorption should remain unchanged and hence these cal- 
ibrations may systematically over-estimate the extinc- 
tion. 

3. PROGENITOR OBSERVATIONS AND DATA ANALYSIS 
3.1. Archival data 

We have searched both the Hubble and Spitzer Space 
Telescope archives, along with the publicly available 
archives of all the major ground based observing facil- 
ities, for pre-explosion images covering the site of SN 
2012aw. A log of all data used is given in Table [TJ The 
deepest, and highest resolution images of M95 found were 
from the Wide-Field and Planetary Camera 2 (WFPC2) 
onboard the Hubble Space Telescope (HST). The site of 
SN 2012aw was observed on multiple occasions between 
1994 Nov and 1995 Jan with the i^439Pv", F555W and 
F8UW filters, and again in the F555W filter in 1995 
Dec, as part of the HST Key Project on the Extragalac- 
tic Distance Scale. HST +WFPC2 observed M95 again 
in F336W and F658A on 2009 Jan 18. M95 was also 
observed with the Subaru Telescope+ SuprimeCAM in 
R on 1999 Jan 28, but as the seeing was > 2.5" in these 
images, they were of no use for constraining the progeni- 
tor. NTT+SOFI (K s ) and VLT+ISAAC (J s ) images of 
M95 were obtained from the ESO archive. A subset of 
the data used is shown in Fig. [2] 

All HST data were downloaded from the Multimission 
Archive at STScfEL The individual WFPC2 F8UW im- 
ages we re combined using the drizzle algorithm within 
iRArQ ((Fruchter fc Hookl 12002ft . although as the offsets 
between the individual exposures were not in an optimal 
subpixel dither pattern, we did not resample the image 
to a smaller pixel scale. Small image shifts were removed 
using the tweakshifts task to improve the quality of 
the alignment. 

The individual VLT+ISAAC images were reduced 
(bad-pixel masked, flat-fielded and sky-subtracted) us- 
ing the ISAAC pipeline (version 6.0.6). The individual 
frames from each night where the SN position was in the 
field of view were then aligned with sub-pixel shifts and 
combined within iraf. The image from the second night 
of ISAAC observations (2000 Mar 27) was used for all 
of the subsequent analysis, as the total on source expo- 
sures are deeper, and the SN position is further from the 
edge of the chip. The NTT+SOFI data from 2006 Mar 
24 were not flat-fielded (as appropriate calibrations were 
not available), but were sky subtracted with the SOFI 

10 http://archive.stsci.edu 

11 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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pipeline (version 1.5.4) before the individual on-source 
frames were combined. A total of 600s was obtained on 
source, with a FWHM of 0.6". For the SOFI data from 
2002 Mar 25, the SOFI pipeline was unable to reduce the 
data, and so the I RAF XDIMSUM package was used for the 
reduction. 

The Spitzcr Space Telescope + IRAC observed M9 5 
as part of the SINGS survey (jKennicutt et al.l [2003), 
however the spatial resolution of these data (1.7" PSF 
FWHM in the shortest wavelength channel) is ill suited 
for progenitor identification. We downloaded the IRAC 
data from the Spitzer Heritage Archive, but could not 
identify a point source at the SN coordinates. We have 
not considered the IRAC data any further in the follow- 
ing. 

3.2. Progenitor identification 

A deep, high resolution i-T-band image of SN 2012aw 
was obtained on 2012 Mar 31 with Gemini+NIRI (Pro- 
gram: GN-2012A-Q-38). To match the resolution of 
the pre-explosion WFPC2 images, we used the ALTAIR 
adaptive optics system to correct for the seeing; as the 
SN was sufficiently bright (mag~13) at the time of the 
observations it was used as a natural guide star for AL- 
TAIR. The fl4 camera was used, which has a 0.05" pixel 
scale over a 51"x51" field of view. As the region around 
SN 2012aw is not crowded, we used on source dithers of 
a few arcseconds, and median combined these to create 
the sky frame for each exposure. The data were reduced 
with the iraf Gemini package, yielding a final combined 
image with a PSF FWHM of 0.2". 

40 sources were identified common to both the NIRI 
and the drizzled WFPC2 F81AW images. The positions 
of these reference stars were measured in both images 
with IRAF PHOT. After rejecting outliers from the fit, 
which either had low signal to noise or centering errors, 
we used 22 sources to derive a geometric transformation 
between the pixel coordinate systems of the WFPC2 and 
NIRI images. A general fit allowing for scaling, a shift in 
x and y, rotation and a skew term was used, with a root 
mean square error in the fit of 53 mas. The pixel coor- 
dinates of SN 2012aw were measured in the NIRI image, 
and transformed to the pixel coordinates of the WFPC2 
image. The transformed positio n of the SN coincided 
wit h the same source id entified bv lElias-Rosa et al"1l2012l 
and lFraser et al. 1120121 The position of this source was 
measured using the three different centering algorithms 
in PHOT, all of which agreed to within 14 mas. The sep- 
aration between the transformed SN position, and the 
progenitor candidate was 27 mas, which is well within the 
total (transformation + SN position + progenitor posi- 
tion) uncertainty of 55 mas. Hence we formally identify 
the source as the progenitor candidate for SN 2012aw. 

3.3. Progenitor photometry 

Photometry was performe d on the WFP C2 images 
with the hstphot package ([Dolphin! [2000a') , which is 
designed specifically for undersampled data from this 
instrument. The pipeline reduced images were first 
masked for bad pixels, hot pixels and cosmic rays us- 
ing the ancillary packages distributed with hstphot. 
After masking, PSF-fitting photometry was performed 
on the individual images, which typically had exposure 



times of 1000 to 2000 s. The SN progenitor candidate 
was detected at a magnitude of F555VF=26.70±0.06 and 
F814iy =23.39±0.02 (in the photometric system as de- 
fined in iDolphml l2()()()'n; . but was not detected in either 
the F33QW or ^4391^ filter images. Based on photom- 
etry of sources detected by hstphot in the F439W im- 
age, we estimate a 5tr limiting magnitude of F439W >26. 
The F814W magnitude for the progenitor, as measured 
in the individual frames, appears to be unchanged within 
the uncertainties over the ~50 day period from the first 
to the last observation. Unfortunately the progenitor is 
too faint to check for variability in the F555W images. 
The i^658iV (Ha) WFPC2 image was also examined, but 
there was no sign of emission at the SN position. 

The SOFI K s and ISAAC Js images were aligned to 
an i ;l 814W / -filter mosaic of the four WFPC2 chips, and 
a counterpart to the optical progenitor candidate was 
found in both images. PSF-fitting photometry was per- 
formed on the ISAAC J s and SOFI K s image with the 
SNOOPY package in iraf. A magnitude was found for the 
progenitor of J = 21.1±0.2 (setting the photometric ze- 
ropoint from the J magnitude of two 2MASS sources in 
the field, which is the dominant source of error). As the 
bandpasses of the J s and J filters are different, we calcu- 
lated colour terms from synthetic photometry of MARCS 
model spectra between 3000 and 4000 K. Over this tem- 
perature range the colour term is < 0.02 mag, which is 
significantly less than our photometric error. For the 
SOFI K s image from 2006, we measured a magnitude of 
K = 19.3 ± 0.4 with PSF-fitting photometry. The un- 
certainty in the measurement consists of 0.36 mag from 
the progenitor photometry, and 0.21 uncertainty in the 
zero point, which was set from seven 2MASS sources in 
the field. As a check of the PSF-fitting measurement, 
aperture photometry was performed on the same image 
with iraf phot, which returned a magnitude of K — 
19.1±0.4. Aperture photometry of the SOFI image from 
2002 yielded a magnitude for the progenitor of 18.9±0.3 
mag. We hence adopt an average value of if =19.1 from 
the two epochs for the progenitor magnitude, with a con- 
servative error of ± 0.4 mag. All NIR m agnitudes are in 
the 2MASS system ([Cohen et al.ll2003| ). 

4. RESULTS AND DISCUSSION 

We used our own Bayesian SED-fitting code (based 
on the Bayesian Infer ence-X nested sampling framework; 
Maund 2012 in prep; lSkilling|[20Cil to compare the ob- 
served progenitor magnitu des to synthetic photo metry 
of MARCS model spectra (|Gustafsson et al.l l2008). The 
models used were for 15 M stars with log g = dex and 
solar metallicity, with temperatures between 3300 and 
4400 K. The SEDs were rebinned by a facto r of 10 prio r 
to the calculation of synthetic photometry ([Pled [20081. 
We fit the progenitor VIJK magnitudes with flat pri- 
ors on T c ff and E(B-V), which was allowed to vary be- 
tween —0.5 < E(B-V) < 2. A match was found with the 
observed progenitor magnitudes along a banana-shaped 
region as shown in Fig. [3] 

The luminosity is ill constrained due to the uncertain 
extinction, and indeed varies within the la contours be- 
tween 5.0 dex at 3550 K to 5.6 dex for a progenitor with 
T eff =4450 K. The SED fit is poorer in K than the J or 
WFPC2 data, as can be seen in Fig. [3] Besides the larger 
photometric errors in K data, the models are also quite 
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sensitive to metalliciity at cooler temperatures, due to 
the TiO absorption which is present in the NIR spectra 
of cool stars. Determining the luminosity of the progen- 
itor from the fC-band magnitude only (as Ak ~0.1 Ay, 
and the bolometric correction to X-band only changes by 
1.1 mag between 3300 and 4500 K) gives a value between 
4.9 and 5.5 dex. We also attempted to fit the progeni- 
tor SED with MARCS spectra of twice solar metallicity 
models, but found the fit to be poorer than for the solar 
metallicity models. 

In Fig. 02 we have also plot ted stellar evolution - 
ary tracks from the STARS code (jEldridge et al.ll2008f l. 
Comparing the luminosity of the progenitor to the end- 
points of the tracks implies a ZAMS mass of between 14 
and 26 M Q . If we consider the hotter progenitors in Fig. 
[U the discrepancy in temperature with the end points 
of the STARS evolutionary tracks becomes more appar- 
ent. However, increased mass loss, either in a binary or 
through rotation, would serve to bring th e end points o f 
these tracks over to hotter temperatures (|Georgvi 2012). 

As the radius of the progenitor can be expressed as 

R=^J L/T^j:, it is easy to calculate the expected radius 

for each point in the HR diagram. We find that even 
for a temperature of 4500 K, and a luminosity of 5 dex, 
the radius of the progenitor is still >500 R©. This is 
sufficiently large that we may reasonably ex pect the pro- 
genitor to give rise to a Type IIP SN (e.g.. |PopovlH993t 
iKasen fc Wooslevl 12009). and so cannot be used to fur- 
ther restrict the region of the HRD where the progenitor 
lies. 

Within the 68% confidence contours, the solar metal- 
licity SEDs favour an extinction which is greater than 
E(B-V) > 0.8 mag, although within 95% (2a), there is a 
solution with a correspondingly poorer fit for E(B — V) = 
0.4 mag. The implication of this is that a significant 
amount of dust could have bee n destroyed in the initial 
phases of the SN explosion (eg. iWaxman fc Draiiiell200Ct 
iDwek et"aT1l2008D . and that the extinction towards a SN 
can not be taken as a proxy for the extinction towards 
the progenitor. We stress, however, that regardless of 
extinction, SN 2012aw potentially arises from one of the 
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highest mass progenitors found t o date. 

As was recently suggested by iWalmswell fc Eld ridge 
(|201lD . significant amounts of circumstellar dust around 
SN progenitors is an appealing solution to the lack 
of high mass red s upergiant progenitors identified by 
ISmartt et all (|2009D . We caution however, that SN 
2012aw is the reddest SN progenitor found thus far, 
and appears to suffer from extinction that is comparable 
to that of the most l uminous Galactic red supergiants 
(|Levesaue et al.l [20051) . but higher than is typical for a 
Type IIP progenitor (from the limited sample with colour 
information). Hence it remains unclear whether circum- 
stellar dust truly is the panacea for the "red supergiant 
problem" . 
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phase [d] Wavelength (A) Wavelength (A) 

Fig. 1. — Follow-up observations of SN 2012aw. Left panel: pseudo-bolometric (BVRI) lightcurve of SN 2012aw, as compared to other 
Type IIP SNe. Centre panel: comparison of early spectra of the same SNe to SN 2012aw (all at rest wavelength). Right panel: Sequence 
of spectra for SN 2012aw over the first two weeks since discovery. 



SN 2012aw 



7 




Fig. 2. — The progenitor of SN 2012aw as detected in archival pre-explosion images. Images are scaled and oriented as indicated, each 
panel is centred on the progenitor position as indicated with tick marks in the F81AW image. In the J-band image the progenitor is 
partially blended with a nearby source, this has been corrected for with PSF-fitting photometry. 
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Fig. 3. — Left panel: posterior probability distribution for reddening and temperature for solar metallicity MARCS SEDs (Contours are 
shown for 68% and 95% credibility intervals). The inclusion of the i-f-band observation (solid contours) provides a significant constraint 
on the possible values of reddening, over the fits that only utilise the observed VI J fluxes (dotted contours). Centre panel: Best fit SED 
(dark red line) for T ef j=3660 K and E(B-V)=1.12), as compared to the observed progenitor magnitudes. The shaded region indicates the 
allowed range of spectra (corresponding to the 68% region in the left panel) with an upper bound of T c g =4490 K, E(B-V)=1.758 mag, and 
a lower bound of T c ff=3500, E(B-V)=0.74 mag. Right panel: Hertzprung-Russell diagram showing the final probability contours for the 
progenitor of SN 20f2aw, together with the STARS evolutionary tracks (labelled with ZAMS mass in Mq) as discussed in the text. 
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TABLE 1 

Pre-explosion images of the SITE OF SN 2012AW 



Date 



Telescope Instrument Filter Exposure (s) 



2009-01-18 

1994-11-29 

-1994-12-19 

2009-01-18 

1994-11-29 

-1995-12-04 

1994-11-29 

-1995-01-16 

2000-03-26 
2000-03-27 

2002-03-25 
2006-03-24 



HST 



WFPC2 



VLT 



NTT 



ISAAC 
SOFI 



F336W 


4400 


F439W 


5000 


F658N 


1800 


F555W 


34130 


F814W 


9830 


J 


2400 




4080 


Ks 


600 




600 



